function meanerr(time,~,R,~,y0,Y)

load('orthbase');

 exactmean=zeros(1,100);
 time(1)=10^-15;

 for ii=1:100
     exactmean(ii)=quadl(@(s) R(s),0,time(ii));
 end
exactmean=y0(1)+exactmean-1/4*time.^2;
hold on
plot(time,exactmean);
plot(time,Y,'r-');
title('Mean')
legend('exact mean','wce mean')
